Atomistics of vapor-liquid-solid nanowire growth 
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Vapor-liquid-solid (VLS) route and its variants are routinely used for scalable synthesis of semi- 
conducting nanowires yet the fundamental processes remain unknown. Here, we employ atomic-scale 
computations based on model potentials to study the stability and growth of gold-catalyzed sili- 
con nanowires (SiNWs). Equilibrium studies uncover segregation at the solid- like surface of the 
catalyst particle, a liquid AuSi droplet, and a silicon-rich droplet- nanowire interface enveloped by 
heterogeneous truncating facets. Supersaturation of the droplets leads to rapid ID growth on the 
truncating facets and much slower nucleation-controlled 2D growth on the main facet. Surface dif- 
fusion is suppressed and the excess Si flux occurs through the droplet bulk which, together with 
the Si-rich interface and contact line, lowers the nucleation barrier on the main facet. The ensuing 
step flow is modified by Au diffusion away from the step edges. Our study highlights key interfacial 
characteristics for morphological and compositional control of semiconducting nanowire arrays. 



Introduction 

A characteristic feature of the VLS route for nanowire 
synthesis is the presence of a catalyst particle, usually a 
droplet, that mediates the mass transfer from the vapor 
phase precursor to the growing nanowire [THS]. The tech- 
nique was demonstrated successfully in the 1960s in the 
context of crystalline whisker growth [1, and has become 
increasingly relevant of late as it offers direct control over 
nanowire composition and morphology, aspects of critical 
importance in several nanowire-based applications [TnlSj . 
The salient features can be readily seen in the schematic 
illustration in Fig. [l^ (inset). A low melting eutectic 
droplet is supersaturated following breakdown of the pre- 
cursor gas on its surface. Past a critical point, crystal- 
lization at the droplet-nanowire interface results in tip 
growth. The droplet then acts as a conduit for ID crys- 
tal growth with diameters that are set by its dimensions. 

The atomic-scale structure, energetics and dynamics of 
the alloyed particle directly influence the morphology and 
composition of the as-grown nanowires, yet the prevailing 
ideas continue to be based on continuum frameworks [12 , 
[M l [T9H24] . As an example, gold-catalyzed growth of sil- 
icon nanowires is a well-studied system yet the morphol- 



* Electronic address: lmupmanyu@neu.edu| 



ogy and composition of the particle-nanowire interface is 
unknown. Basic crystal growth parameters such as the 
size dependent equilibrium composition Xg- above which 
the droplet rejects the excess silicon onto the nanowire 
remain unknown; continuum models models ignore the 
effect of segregation and related interfacial phenomena. 
Such microscopic features modify the droplet morphology 
through the balance of interfacial tension at the contact 
line which in turn dynamically regulates the chemical po- 
tential difference between the two phases, i.e. the driv- 
ing force for growth. Likewise, the atomic-scale processes 
that sustain the growth of the nanowire remain unknown. 
The steady-state is limited by reduction of the precursor 
gas on the droplet surface [25 and the 2D nucleation- 
limited step flow at the droplet-nanowire interface [27l? 
]. Understanding the interplay with the dynamics that 
feeds the growth of the nuclei is crucial for compositional 
control of the as-grown nanowire I28H3Q] and requires 
an atomistic understanding of the near-equilibrium par- 
ticle behavior. 

We focus exclusively on a small diameter 2R = 10 nm, 
(111) silicon nanowire growing isothermally above the 
eutectic temperature at T = 873 K. Both experiments 
and ab-initio computational techniques are handicapped 
in accessing the spatio-temporal scales associated with 
the energetics and 3D dynamics during growth [31^33]. 
We capture these processes in their full complexity by 
employing atomistic simulations based on an angular 
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FIG. 1: Phase diagram of the model Au-Si system: (a) Bulk phase diagram for the Au-Si binary system predicted by 
the model inter-atomic interactions employed in this study (solid black line) and the size correction to the model diagram 
Xg^{2R) = 0.46 (solid red circle) based on doubly-capped SiNW computations shown in (b). The solid gray line is the 
experimental bulk phase diagram shown for comparison, (inset) Schematic illustrating the VLS route for SiNW nanowire 
synthesis, (b) (left) Atomic configuration of an equilibrated 10 nm long SiNW capped with catalyst particles at both ends, 
(right) Mid-section slice through the droplet-nanowire system that reveals the shape of the droplets and the faceted nature 
of the solid nanowire. Light yellow is solid silicon, and red and bright yellow indicate liquid Si and Au, respectively. For 
consistency the color scheme for the atomic plots and the density profiles is preserved in the remaining figures. 



embedded-atom-method (AEAM) framework for describ- 
ing the hybrid metallic-covalent interactions [34]. The 
computed surfaces of equilibrated isolated droplets are 
ordered and Si-rich, yet decorated by a submonolayer 
of Au. The surface structure leads to sluggish atomic 
mobilities such that most of the excess Si flux during 
growth is directed through the droplet bulk. Studies on 
droplet-capped SiNWs show that the solid-liquid inter- 
face is again Si-rich and the nanowire sidewalls are dec- 
orated by Au monolayers. Of note is the interface mor- 
phology composed of a main facet that truncates into 
smaller staircase-like, stepped facets at the contact line. 
The latter grow rapidly in a ID fashion while crystal- 
lization on the main facet proceeds via nucleation of 2D 
islands which then grow via step flow. Interestingly, the 
motion of the steps is limited by Au diffusion into the 
droplets and/or the nanowire sidewalls. Studies with 
varying levels of droplet supersaturation yield the nucle- 
ation barrier which is lower than that for homogeneously 
nucleation on a planar Si(lll)-AuSi interface, highlight- 
ing the combined role of segregation, atomic diffusivities, 
and interface morphology in regulating nanowire growth. 



Results 

A. Size corrected liquidus composition 

We first quantify the equilibrium state of the alloy sys- 
tem by extracting the liquidus composition, X^-{2R^ T = 
873 K). We capture the effect of the solid-liquid inter- 
face and the contact line by extracting the equilibrium 



state of a 2R = 10 nm diameter, (111) SiNW capped 
by a droplet at each end. The geometry allows the 
nanowire to act as a source/sink as the droplets adjust 
to the size-and morphology-dependent liquidus composi- 
tion. The starting configuration consists of a nanowire 
with a hexagonal cross-section bounded by {112}-like 
sidewall facets (Methods). In particular, we have per- 
formed separate computations for SiNWs with {110} and 
{112} sidewall facets. The choices are motivated by past 
experiments that have shown that the sidewalls can be 
a mixture of {112} and {110} facets [TS', l35]. Figure [i]d 
shows the atomic configuration of the equilibrated SiNW 
system, in this particular case bounded by {110} facets. 
We immediately see the formation of an Au submono- 
layer on the sidewalls, consistent with prior experimen- 
tal observations [36 . The cross-section maintains its six- 
fold symmetry for both {110} or {112} family of facets 
implying that the two surfaces are energetically quite 
close; we do not observe their co- exist ence p^ [35]. This 
is not surprising as the thermodynamic shape is inher- 
ently size dependent [37J and the diameters of the synthe- 
sized SiNWs are much larger. The {112} facet is unsta- 
ble for Au coverages typically observed during nanowire 
growth[13 , yet we do not observe axial serrations due 
to the small nanowire lengths. We also cannot rule out 
the effects of surface stresses at these small sizes [38 as 
they can change the relative stability between {110} or 
{112} family of facets. The average composition within 
the interior of the capping droplets is Xg-{2R) = 0.46 
and is indicated on the phase diagram in Fig. [1^. The 
size-corrected supersaturation with respect to the bulk 
liquidus composition is Ax = Xg-{2R) — Xg- ^ 1% and 
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FIG. 2: Surface properties of equilibrated AuSi droplets: 

(a) Midsection through an atomic plot of a 32,000-atom AuSi sys- 
tem equihbrated at the model liquidus composition X^^ = 0.46. (b 
and c) Radial variation of (b) the ensemble averaged and normal- 
ized density profile p(r) / p for the combined system (solid dark line) 
and for Au and Si, as indicated, and (c) ensemble averaged diffu- 
sivities of Au and Si, ^Ati/Si(^) extracted from equilibrium MD 
simulations of the droplets. The dotted lines in (c) correspond to 
the bulk diffusivities extracted from simulations on an equilibrated 
bulk alloy system (see Methods). 



the equivalent undercooling At ~ 18K. 



B. Surface segregation, ordering and kinetics 

The composition of the capping droplets is non- 
uniform and we decouple the surface effects by first 
studying isolated droplets at Xsi = Xg-{2R). Atomic 
configuration of a midsection through an equilibrated 
droplet is shown in Fig. [2^. The non-uniform compo- 
sition arises mainly due to surface segregation. Details 
are extracted using normalized density profiles {p{r)/p) 
plotted as a function of the distance r from the droplet 
center of mass (Fig. |2]3). The overall profile (dark solid 
line) has a well-defined peak at the surface suggestive of 
a solid-like surface structure. The multiple alternating 
peaks in the Au and Si profiles indicate a Si-rich sub- 
surface decorated by a submonolayer of Au. While the 
Si segregation is consistent with experiments on eutectic 
AuSi thin films [39], the Au submonolayer has not been 
reported earlier. 

We further isolate the effects of surface curvature via 
segregation studies on thin AuSi films and the results 
are summarized in Supplementary Figure SI. At the eu- 
tectic point, Xsi ^ 0.31 and T = 590 K, we see ripples 
in the density profiles that extend to ^7 layers and de- 
cay rapidly into the bulk. The density variations are 
larger than those usually seen on liquid surfaces and sug- 
gest a compositionally ordered pre-frozen surface quali- 
tatively similar to the surface ordering seen in eutectic 
AuSi thin films [39]. The ordering decreases with increas- 
ing Xsi and temperature, and the extracted density pro- 
files are similar to those observed in the liquid droplets. 



i.e. the surface curvature has little effect. The surfaces 
are always decorated by a Au submonolayer. It is ex- 
pectedly less ordered compared to monolayers on crys- 
talline Si surfaces, and this feature likely makes it dif- 
ficult to be resolved by X-ray based surface characteri- 
zation techniques [39]. The discrepancy can certainly be 
an artifact of the EAM framework used to describe the 
metallic bonds [7 ; it underestimates the surface tension 
of liquid Au by 20% and can therefore artificially favor 
the formation of the Au submonolayer. Nonetheless, it 
is natural to expect that the stability of Au monolay- 
ers on crystalline Si surfaces [41 also extends to partially 
crystalline Si surfaces. As confirmation, we still see the 
Au decoration in segregations studies with modified Au- 
Au interactions using a charge gradient approach fit to 
the surface tension of liquid Au[7 (Supplementary Fig- 
ure S2 and Supplementary Methods). The Au coverage 
is slightly reduced in extent yet the qualitative trends 
remain unchanged. The effect of surface structure ex- 
tends to the diffusive kinetics as well. We quantify the in- 
terplay by monitoring the atomic mean-square displace- 
ments (MSD) that yield the atomic diffusivities Dj^^^si 
(Supplementary Figure S3). The radial variation is plot- 
ted in Fig.|2]3 (Methods). For both Au and Si, the surface 
diffusivity is less than half of the bulk value and further 
highlights the dramatic effect of the surface structure. 



C. Equilibrated nanowire-particle system 

The substrate is explicitly included by mating an equi- 
librated droplet- SiNW system onto a dehydrogenated 
Si(lll) surface (Methods). Figure |3^ and Supplemen- 
tary Figure S4 show the atomic configurations of the 
nanowire system with {110} and {112} sidewall facets, 
respectively. A stable Au submonolayer is evident on the 
nanowire sidewalls (right. Fig. [3|i). The surface segrega- 
tion is similar to that on isolated droplets (Fig. [2]3), i.e. 
the substrate has a minor effect. There are differences in 
the structure of the Au submonolayer for the two classes 
of sidewall facets. We see Si segregation at the droplet- 
nanowire interface within liquid layers adjacent to the 
interface and parallel to (111) planes. Ripples in the 
in-plane density profiles within the droplet extracted as 
a function of distance from the interface z extend to a 
thickness of ~ 5 A (Fig. [sj)). The sharp Si and Au peaks 
adjacent to the interface at 2; ~ 1.5 A (arrow) indicate 
a Si-rich layer with some in-plane order. The remain- 
der of the segregated layer is mostly liquid characterized 
by a smaller and broader Si-peak at z ~ 2.5 A (arrow), 
followed by a Au-rich region (z ~ 4A) which ultimately 
settles into the bulk composition. The compositional dis- 
tribution for the entire nanowire-droplet system is sum- 
marized schematically in Fig. [Sj). 

Surface plot of the solid nanowire sans the the liquid 
atoms yields insight into the interface morphology (mid- 
dle panel. Fig. [^). Clearly, it is not flat. The (111) 
facet is dominant yet it is truncated along the periph- 
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FIG. 3: Equilibrated droplet-nanowire-substrate system: (left) Equilibrium atomic configuration of a 2R = 10 nm diameter (111) 
silicon nanowire with {110} sidewall facets. The wire is grown on a (111) silicon substrate and is capped by an AuSi catalyst particle at 
the liquidus composition, X^^ = 0.46. (middle, inset) The surface plot viewed along the [111] direction showing the morphology of the 
solid SiNW that abuts the solid-liquid interface. The probe radius is fixed at 1.2 A. Schematic illustration of one of the corners (arrow). 
The disordered region is shaded gray, (right) Atomic-plots of the midsections along the indicated directions showing the distribution of 
the {113} and {120} family of truncating facets fianking the interface, (b) The in-plane density profile py averaged over each liquid bilayer 
abutting the <111> interface [shaded strip in (a)] as a function of the distance z. The arrows indicate peaks consistent with Si segregation 
within the bilayer. (c) Schematic illustrating the interplay between the segregation profile and the interfacial energetics that maintains the 
local equilibrium at the contact line and the droplet shape, mainly angles a and /3, main facet length A(m), and truncating facet lengths 

^{113} ^{120}- 



ery next to the contact line. These truncated regions are 
2-4 atomic diameters wide and their length varies non- 
uniformly along the contact line. Multiple midsections 
reveal a stepped structure yet they are not smoothly 
curved; rather they correspond to two distinct orienta- 
tions inclined at /3 ^ 32° and /3 ^ 50° to the (111) side- 
wall facet (vertical), which we identify to be the {113} 
and {120} family of planes. Supplementary Figure S4 re- 
veals a similar morphology for SiNWs with (112) sidewall 
facets with the exception that the truncating facets con- 



sist of {113} and {111} family of planes. The six facets 
alternate between these two inclinations along the con- 
tact line, as indicated in the surface plot and the midsec- 
tions (Fig.[3|i, right). The facets are also observed in the 
doubly-capped nanowire system shown in Fig. ^p. The 
six edge corners where truncating facets intersect are in- 
herently rough as they have to absorb the differences in 
the facet inclinations, as illustrated schematically. The 
roughening can also be seen in the surface plot in that the 
step morphology kinks at the corner edges. The droplet 
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shape is therefore asymmetric since the angles a and ^ 
that fohow from the Young's balance along the contact 
line are shaped by the facet energetics and geometry. 
The asymmetry is also evident in the two midsections: 
a ^ 11° at the {113} truncating facet and changes to 
a ^ 15° at the {120} truncating facet. 

The presence of nanometer-wide stable truncating 
facets is consistent with recent reports in several 
semiconducting nanowire systems, including the Au-Si 
svstem[42l-[44]. Limitations of the model system notwith- 
standing, such a complex 3D morphology composed of 
multiple families of truncating facets is similar to that 
recently reported in [0001] sapphire nanowires[42 , sug- 
gesting that this is perhaps a general feature during VLS 
growth. For the specific case of Si, it is interesting 
that the {113} and {111} facets routinely decorate the 
sidewalls during the complex sawtooth faceting in much 
larger diameter (111) SiNWs [I3l [35 , suggesting that 
they perhaps nucleate first at the solid-liquid interface. 



D. Nanowire growth kinetics 

The heterogeneous nature of interface facets has im- 
portant ramifications for nanowire growth which we now 
study computationally. We simplify the catalysis by 
directly abstracting individual silicon atoms onto the 
droplet surface to a prescribed level of supersaturation, 
and employ classical molecular dynamics to study the 
non-equilibrium response of the entire droplet-SiNW- 
substrate system (Supplementary Methods). Figure [I] 
shows the structural evolution at and around the inter- 
face as well as the in-plane density profile p\\{z) for a 
droplet with initial composition Xsi = 0.48 (A^ ~ 36 K). 
The initial growth stage for t < 0.4 ns is marked by layer- 
by-layer crystallization on the truncating facets. The cir- 
cumferential growth is likely aided by their small widths 
and the proximity to the Si-rich contact line (circled. 
Fig. [4^). We see ID lengthening of kinks associated 
with the stepped structure of these facets. Interestingly, 
the growth occurs preferentially away from the corner 
edges implying that they favor dissolution. For details, 
see Supplementary Movie 1. The extent of truncation 
changes due to their growth, consistent with the oscilla- 
tory growth wherein the facet widths decrease between 
bilayer additions on the main facet [42- 44 . We also see 
smaller clusters of 2-3 Si atoms that continually form and 
dissolve on the main facet yet the density profile normal 
to the interface remains unchanged (bottom row), i.e. the 
crystallization is limited to the truncating facets. 

The t = 0.8 ns configuration plotted in Fig. ^ shows 
that, barring the corner edges, the growth has advanced 
by almost a full layer on the truncating facets. The 
droplet is still supersaturated and it now drives the nu- 
cleation on the main facet. It is worth pointing out that 
structurally the growth on the truncating facets and on 
the main facet is essentially independent, since a 2D is- 
land reaching the edge of one facet has no nearest neigh- 



bors with an island on the adjacent facet. They require 
the influence of the much weaker next nearest-neighbor 
interactions to extend across. We observe early stages 
of classical 2D layer-by- layer growth wherein several nu- 
clei form and attempt to overcome the nucleation bar- 
rier. Observe that the nuclei form slightly away from the 
center of the main facet and their location is correlated 
with the corner edges. The disorder together with the 
enhanced Si-segregation at the contact line likely aids 
the nucleation on the main facet. The diffusion along 
the surface is indeed slower yet the bulk liquid can easily 
provide a pathway. The precursor cluster that eventually 
becomes supercritical is circled in Fig. [4]3. Figures [4]3-|4^ 
reveal the key features of the subsequent bilayer growth 
on the main facet, i.e. it remains layer-by-layer and dom- 
inated by a single growing nucleus. 

The density profiles evolve accordingly. The intensity 
of the two silicon peaks increases at the expense of the Au 
peak. A sharp peak in the adjacent liquid layer indicates 
build-up of Si while the third layer becomes Au-rich (ar- 
row in Fig.[4]i). Evidently, Au diffuses rapidly away from 
the interface. At t = 2 ns, the first layer has not com- 
pletely formed yet the step flow has slowed down due to 
the low supersaturation left in the droplet, Xsi = 0.463. 
Note that there is negligible growth at the edge corners 
(circled. Fig. [4^). We see Au and liquid Si within the 
mostly crystallized Si layer although the liquid Si atoms 
exhibit in-plane order. Eventually, some of the Au atoms 
diffuse away into the liquid droplet and aid further crys- 
tallization (not shown). 

We quantify these observations by monitoring the crys- 
tallized volume at the interface V{t) and the droplet com- 
position ^^^(t) (Fig.|4^). Crystallization kinetics on the 
truncating and main facets are plotted separately. Ini- 
tially for t < 1 ns the increase in V{t) is entirely due 
to crystallization on the truncating facets. A linear fit 
yields a growth rate of ^ 10.5 cm/s. The growth is 
instantaneous indicative of a small nucleation barrier. 
While the barrier can be artificially reduced due to the 
larger supersaturations in the computations, note that 
the Si-rich contact line is an easy source for Si adatoms. 
Also, as in the experiments, the facet widths rarely ex- 
ceed a few nanometers such that the critical nuclei sizes 
are likely larger than the facet widths [45 . The circum- 
ferential ID growth is also aided by the stepped structure 
of these facets and we quantify this effect by extracting 
Si adatom enthalpies Had on the {111} and the {113} 
family planes (Methods). Not surprisingly, the close- 
packed {111} facets have a significantly higher enthalpy 
than the stepped {113} facets, AHad ^ 2.0 eV/adatom. 
The adatom diffusivity on the {113} facets is the sup- 
pressed since the atom makes fewer diffusion hops be- 
fore another atom crystallizes in close proximity. The 
adatoms can readily form clusters and therefore lower 
the nucleation barrier, which in turn can induce rapid 
ID growth. The role of such structural effects in intro- 
ducing kinetic anisotropics during faceted crystal growth 
is well-established [46^. 
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FIG. 4: Nucleation and growth kinetics at the nanowire-droplet interface: (a-e) Atomic configurations in the vicinity of the 
interface subject to a supersaturated droplet with initial composition Xsi = 0.48. (top row) Liquid bilayer abutting the droplet-SiNW 
interface, and (middle row) top and side views of the solid nanowire. The light and dark blue colors indicate crystallized Si atoms on the 
truncating and main facets, respectively. The instantaneous supersaturation is indicated on the top of each panel. The circled region in 
(a) highlights an instance of preferential growth at the center of the truncating facet while the circled region in (b) indicates precursor 
fluctuations that lead to the formation of a critical nucleus on the main facet. Details can be seen in Supplementary Movie 1. The bottom 
of each panel shows in-plane density proflle p\\(z). (f) Schematic illustration of the evolution of the interface before and during growth 
via step flow. See text for discussion, (g) Temporal evolution of crystallized net Si volume (red circles) composed of crystallization on 
the edge and main facet (light and dark blue circles, respectively), as well as the droplet composition (solid green squares). The linear flt 
(solid red line) used to extract the layer growth rate due to post-nucleation step flow on the main facet (h) Size of the nuclei observed on 
the main facet. The red curve is the evolution of the critical nucleus as indicated, while the other curves correspond to evolution of the 
smaller peripheral nuclei. 
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At t ^ 1.0 ns, the growth rate on the truncating 
facets decreases, either due to site-saturation or the on- 
set of the oscihatory growth mode. The growth on the 
main facet, non-existent so far, dramaticahy increases. 
This is indicative of the growth of a supercritical nu- 
clei and is again consistent with our qualitative observa- 
tions (Fig. [4j3-|4^). In order to ascertain if the growth is 
layer-by- layer or kinetically rough, we monitor the size 
evolution of the visibly identifiable nuclei on the main 
facet (Fig.[4[i). We see a dramatic increase in the size of 
one of the nuclei while the other nuclei fiuctuate about 
a relatively smaller size (< 35 atoms). At longer times 
t ^ 10 ns, one of these smaller nuclei shrinks while the 
other two persist without any appreciable increase in size. 
Clearly, one 2D island dominates the layer growth and 
the encompassing steps are the main contributions to 
surface roughening. Moreover, the step orientations re- 
fiect the symmetry of the growth direction while that is 
not the case for the fiuctuating subcritical islands. The 
size evolution shows that even though the driving force 
is larger than that in experiments, kinetic roughening 
or the increase in step density resulting from growth on 
the main facet is non-existent [47 . The persistence of 
the subcritical nuclei can signify that the supersaturation 
is on the borderline between single nucleus and polynu- 
clear growth; then at lower supersaturations we expect 
the growth of the nanowires to be universally by the sin- 
gle nucleus mechanism. Following nucleation, the growth 
is steady-state as the volume increases and composition 
decreases. A linear fit in this regime yields the over- 
all growth velocity, v ~ 6.3cm/s. The growth rate also 
includes contributions from the truncating facets. Sub- 
tracting the former from the overall growth rate yields 
the step velocity on the main facet, Vg ~ 3cm/s. Fig- 
ure |4]F is a schematic summary of the morphological fea- 
tures implicated in the overall growth. 



E. Effect of supersaturation 

In order to make contact with experimental scales, 
we quantify the nucleation barrier on the main facet 
by performing computations with decreasing initial su- 
persaturation. Representative results for Xsi = 0.463 
(Ax ^ 5.4 K) are shown in Fig. [sj The truncating facets 
again grow rapidly and circumferentially (Fig. [5|i). The 
growth is instantaneous yet slower, and the growth on the 
main facet is non-existent (Fig. jsj)). The crystallization 
is again preferentially at the centers of the truncating 
facets, presumably due to a ID interplay between crys- 
tallization at the facet center and melting at the rough 
corner edges (Supplementary Movie 2). At t ^ 3.0 ns, we 
see transient crystallization of subcritical nuclei on the 
main facet yet there is no supersaturation left for their 
growth. Then, Xsi = 0.463 represents a critical initial 
composition below which there is no nucleation on the 
main facet, i.e. nanowire growth can only occur at larger 
supersaturations. 



Increasing the initial composition in the range Xsi > 
0.485 (At > 45K) decreases the dead-zone that precedes 
the nucleation on the main facet. Multiple stable islands 
contribute to the growth of a single layer such that the 
step density increases and there is a gradual transition to 
a kinetically roughened regime wherein the second bilayer 
begins to crystallize before the complete growth of the 
first bilayer. This is evident in configurations for Xsi = 
0.485 shown in Supplementary Figure S5. Several Au 
atoms remain confined within the first layer. They likely 
retard the step fiow and that triggers the onset of a 3D 
growth mode. Midsections reveal additional Si diffusion 
down the sidewalls that further aids the absorption of the 
supersaturation (inset). 

Figure [6^ shows the driving force dependence of the 
overall growth velocities for varying supersaturations. 
The velocities are linear fits to the initial crystalliza- 
tion rates. In the small driving force limit Xsi < 0.48 
(At < 36 K), the velocity increases rapidly. Results of 
similar computations on a planar Si(lll)-AuSi interface 
are also plotted for comparison (see Methods). Crys- 
tallization at the planar interface entails homogeneous 
nucleation and is always slower at comparable driving 
forces, suggesting a higher nucleation barrier. Separate 
plots for the contributions of the truncating and main 
facets are plotted in Fig. [BfD. The velocity of the trun- 
cating facets increases linearly with supersaturation and 
indicates that its growth is not nucleat ion-controlled, at 
least for undercoolings as low as At = 5.4 K. The av- 
erage kinetic coefficient is ^ 3.0mm/(sK) and is of the 
order of kinetic coefficients in binary alloys [48]. 

The crystallization rate on the main facet increases 
non-linearly with the supersaturation signifying classi- 
cal nucleation-controlled growth. It is possible that at 
much smaller driving forces the growth is linear as the fi- 
nite facet size requires equilibrium between melting and 
nucleation [32]. although it is not clear if this applies to a 
fully faceted interface. Accordingly, we fit the variation 
to a functional form v = 'Uq exp(— AG/Zc^T), where the 
exponential factor is the probability that a spontaneous 
fiuctuation will result in a critical nucleus, and vq is the 
2D growth rate of the supercritical nucleus that is almost 
independent of the supersaturation. The nucleation bar- 
rier varies as AG oc 1/Ax and the exponential fit to 
the extracted velocities yields a value AG = 4meV/Ax 
at T = 873 K. A similar fit to the planar interface data 
yields a barrier that is 63% larger and is consistent with 
estimates based on nucleation of a 2D island (Supple- 
mentary Discussion 1). At larger driving forces charac- 
terized by kinetically rough 3D growth, the velocity in- 
creases almost linearly with supersaturation. The plot 
also shows results of computations on undersaturated 
droplets {Xsi < X^-). At small driving forces, the nu- 
cleation and growth dynamics associated with the melt- 
ing kinetics is completely reversed; the truncating facets 
melt first as the energy of formation of a vacancy on 
these facets is lower, and eventually the main facet melts 
preferentially from the corner edges. The layer- by-layer 
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FIG. 5: Growth kinetics at small supersaturations: (top row) Atomic configurations of the solid nanowire growth front subject 
to an initial droplet supersaturation of Xsi = 0.463. The elapsed time and the current supersaturation are also indicated. For details, 
see Supplementary Movie 2. (bottom) The temporal evolution of the crystallized volume V(t) and the composition within the droplet, 
^Si(t)- (inset) Evolution of the crystallized volume on the main and truncating facets. 



melting velocities are smaller compared to those during 
growth due to the different morphological features impli- 
cated in the two processes, which by itself is important as 
the growth is a steady-state between growth and melting 
kinetics [32^ .45^. 



for both short-range atomic hops into adjacent liquid bi- 
layer and long-range diffusion towards the droplet sur- 
face. The diffusion is markedly along the bulk and also 
towards the nanowire sidewalls; the latter represents a 
direct link between nanowire growth and Au decoration 
on the sidewalls. 



F. Mechanistic insight during growth 

Maps of atomic displacements over prescribed time 
intervals corroborate the observed trends. Figure [7^-d 
shows a t = 0.5 — 1.0 ns map for select Si and Au atoms 
within the Xsi = 0.48 computation. The displacements 
correspond to newly crystallized Si atoms (Si-map) and 
Au atoms initially at the interface (Au-map). The 2D 
step flow results in linear decrease of the droplet compo- 
sition from Xsi = 0.474 to Xsi = 0.469 (Fig. [4^). The 
Si crystallization is mostly from adjacent bilayers and 
driven by displacement chains. Surprisingly, we also see 
long-range fast diffusion from as far as the droplet sur- 
face (Fig.[7^). The flux is mostly through the amorphous 
bulk of the droplet and we see negligible surface diffusion 
(Fig.jrjD). The step flow is aided by Au diffusion into the 
droplet (Figs. [tJ^) and also along the interface directed 
away from the moving steps (Fig. [7]i). We see evidence 



Discussion 

The atomic-scale structure and dynamics during VLS 
growth of SiNWs is clearly sensitive to segregation at 
the surface and the droplet-nanowire interface. The ro- 
bust surface layer is expected to be stable under reac- 
tor conditions where the droplet is exposed to the pre- 
cursor gas atmosphere. The Au submonolayer together 
with the inhibited surface diffusion renders the droplet 
surface in immediate contact with the vapor sufficiently 
Au-rich to drive the precursor breakdown over large 
pressure and temperature ranges. The heterogeneously 
faceted interface morphology impacts all aspects of the 
growth process [43^ and requires reconsideration of cur- 
rent growth models. The truncating facets induce varia- 
tions in the Si segregation and droplet morphology along 
the contact line. The relative disordered corner edges im- 
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FIG. 6: Supersaturation dependence of nanowire growth rate: (a) The extracted linear growth velocity v[= 1/{7tR^) dV/dt] of 
the nanowire as a function of the initial supersaturation, Ax = ^Si ~ ^si' equivalent undercooling based on the phase diagram in 
Fig. [l^ is also indicated. For comparison, growth velocity of a planar Si(lll)-AuSi solid-liquid interface is also plotted (dotted line) as 
a function of initial supersaturation Ax(R = oo) within the AuSi phase. The gray shaded areas represent kinetically roughened growth 
regime induced by large supersaturations. (b) Expanded view of the small supersaturation growth data corresponding to the boxed region 
in (a). The contributions of the two classes of facets and the relevant curve fits are shown as separate plots. 



pact interface nucleation and growth. They hkely serve 
as efficient sources of steps for rapid crystallization on 
the truncating facets, analogous to the screw dislocation- 
based spiral crystal growth [49^. At much lower supersat- 
urations, the dynamics may involve a competition be- 
tween nucleation at the facet centers and melting at the 
edge corners. Nevertheless, the barrier is clearly much 
smaller than that on the main facet. The latter is clearly 
the rate limiting event. The growth of the truncating 
facets also decreases the extent of these ordered regions, 
and they can facilitate mass transfer from the Si-rich con- 
tact line onto the main facet and thereby lower the barrier 
compared to an equivalently supersaturated planar Si- 
AuSi interface; the computations show that this is indeed 
the case. The extracted nucleation barrier AG allows us 
to establish nucleation- and catalysis-controlled regimes. 
Under typical growth conditions, the atomic incorpora- 
tion rate is of the order of a single Si atom per millisecond 
(Supplementary Discussion 2). Equating this incorpora- 
tion rate to the nucleation controlled growth rate, the 
extracted barrier AG = 4meV/Ax yields a supersatu- 
ration of ^ 0.002 for isothermal nanowire growth at 
T = 873 K, or a critical undercooling of A^ ^ 3.6 K. The 
growth is then nucleation-controlled for At < A^ and 
catalysis-controlled otherwise. Solidification from melts 
typically occurs at undercooling of less than a Kelvin [50 , 
and if that is the case during VLS growth of nanowires, 
it must nucleation-controlled. We should emphasize that 
since the barrier is expected to be nanowire size depen- 
dent, the interplay between nucleation and catalysis can 



be quantitatively different at larger sizes and is a focus 
of future studies. 

An interesting coda is the transition to kinetically 
rough growth mode at larger supersaturations is intrigu- 
ing as the step flow is increasingly limited by Au-diffusion 
away from the step-edges. The inter-island Au atoms 
are stable over the MD time scales (< 10 ns), and al- 
though we cannot rule out their slow diffusion out of the 
nanowire eventually, there is growing evidence of catalyst 
particle trapping within the as-grown nanowires [5 1. De- 
tailed understanding of this regime opens up the possibil- 
ity of doping nanowires to non-equilibrium compositions 
during their synthesis [16^, a scenario worthy of detailed 
explorations by itself. 



I. METHODS 

Model system: We use classical inter-atomic poten- 
tials to model the eutectic AuSi system. Hydrogena- 
tion effects are ignored due to hydrogen desorption and 
Au-surface passivation at these elevated temperatures. 
The potentials can be a limitation as they represent ap- 
proximations to the electronic degrees of freedom. Al- 
though the AEAM model system is fit to several exper- 
imental and first-principles data and is based on estab- 
lished frameworks for both pure Au and Si, its ability to 
make quantitative predictions, notably interface proper- 
ties, may be limited. For example, the pure Au potential 
underestimates the liquid free energies^ while the pure 
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FIG. 7: Diffusion pathway during nanowire growth: (a-d) Cumulative displacement maps for select (a-b) Si and (c-d) Au atoms 
over a 0.5 ns time interval following nucleation, for initial composition Xsi = 0.48. Both side (a, c) and top (b, d) are shown. The Si 
displacements correspond to crystallized atoms within the time interval. The Au displacements are for atoms initially at the solid-liquid 
interface. Within each map, the arrow lengths are proportional to the net displacement. Atomic configurations of midsections (side views) 
and bilayers (top views) are superposed within the maps to facilitate visual location of the crystallizing Si atoms (shaded blue). 



Si potential is limited in describing surface reconstruc- 
tions on crystal facets [52 . Nevertheless, validation stud- 
ies presented here show that the predictions of the model 
potential are in agreement with several alloy properties. 
The AEAM framework reproduces the primary feature 
of the binary phase diagram, i.e. a low melting eutec- 
tic. The noteworthy deviations are an underestimation of 
the eutectic temperature Te (590 K compared to 636 K) 
and overestimation of the eutectic composition X^- (31% 
compared to 19%). We have also performed contact an- 
gle studies of eutectic and hypereutectic AuSi droplets on 
Si(lll) surface using semi-grandcanonical Monte-Carlo 
(SGMC) simulations (Supplementary Methods). The fi- 
nal configuration of a 2R = 10 nm size droplet on a 
Si(lll) surface at the growth temperature is shown in 
Supplementary Figure S6. The contact angle = 143° is 
close to the experimental value [1, = 136°. 

Equilibrium liquidus concentration for 2R = 
10 nm SiNW: The initial structure consists of a pristine 
2R = 10 nm (111) SiNW with {110} or {112} family of 
lateral facets capped by two AuSi particles. The Si(lll)- 



AuSi interfaces are initially fiat. The SiNW length is 
^ 8 nm, large enough such that the interfaces do not in- 
teract. The atomic structure of the particles is generated 
randomly at a density corresponding to liquid AuSi and 
the initial composition is hypoeutectic, Xsi = 0.25. The 
entire system consists of ~ 64000 atoms and is relaxed 
using SGMC computations. A local order parameter is 
employed to demarcate the liquid and solid Si atoms (54], 
The order parameter is an extension of the pure Si case 
based on the local 3D structure around each Si atom. 
Additional validation of the surface segregation is done 
using charge-gradient corrections to the Au-Au interac- 
tions that correctly reproduces the liquid Au surface ten- 
sion (Supplementary Methods). 

Equilibrium studies on isolated AuSi particles 
and thin films: Density profiles of the equilibrated 
droplets are averaged over bilayer thick spherical shells 
around the center of mass of the droplet, and the ensem- 
ble average is over 1000 different atomic configurations. 
In order to quantify the effect of droplet surface curva- 
ture, the density profiles are also extracted for thin films. 
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We employ a thin film of thickness H = 10 nm, free at 
both surfaces and periodic in the film plane. SGMC 
simulations are used to equilibrate the initially liquid 
film. Supplementary Figure SI shows the equilibrated 
configuration for a 27000-atom simulation and the cor- 
responding normalized density profiles as a function of a 
distance from the surface h — H extracted at the bulk 
eutectic composition and temperature, Xsi = 0.31 and 
Te = 590 K. Layering of both Au and Si is clearly vis- 
ible and a Au monolayer decorates the Si-rich surface. 
Supplementary Figure SI also shows the density profiles 
for thin films at the growth temperature selected for this 
study (T = 873 K) and for size-corrected liquidus point, 
(T = 873 K, Xsi = 0.46). Comparison with the profile at 
the eutectic point clearly shows a decrease in the thick- 
ness of the solid surface layer with increasing temperature 
and Si supersaturation. 

Au and Si difFusivities within the droplet: MD 
simulations are performed on the equilibrated droplets 
to extract the ensemble-averaged MSD of both Au and 
Si (Supplementary Methods). A net simulation time of 
a few tens of picoseconds is sufficient to extract statis- 
tically meaningful data. For each atomic species at an 
initial distance r from the mass center of the particle, the 
MSD increases linearly with time and the slope yields the 
diffusivityl)(r). For atoms in the bulk and on the surface, 
the slope is different. The MSDs for Au and Si atoms ini- 
tially at the surface of the 2R = 10 nm Au54Si46 droplet 
are plotted in Supplementary Figure S3. The computa- 
tions on bulk eutectic AuSi alloy yield a diffusivity of 
Dau = 1.54 X 10-^mVs at 873 K; the value is in excel- 
lent agreement with past experiments on Au diffusivity 
in eutectic AuSi alloys [55 , Dau ^ 1-7 x 10~^m^/s. 

Equihbrium state of the droplet-nanowire- 
substrate system: The initial configuration is gener- 
ated by taking one half of the double droplet configu- 
ration and placing it on a (16.0 nm x 15.4 nm) Si(lll) 
substrate. The dimensions of the substrate are the same 
as in the contact angle studies. The bottom bilayer is 
fixed. The exposed area is decorated by a Au monolayer 
with a structure similar to the classical x ^/3R30° 
monolayer structure p , and the entire system is equili- 
brated using canonical Monte-Carlo simulations followed 
by MD simulations. 

Nanowire growth kinetics: The catalytic incorpo- 
ration of individual Si atoms takes place over much larger 
time-scales, and we expect the interface growth to be 
orders of magnitude faster. Therefore, there is no addi- 
tional Si flux imposed on the droplet during the course of 
the MD computations and we investigate the effect of an 
initially prescribed driving force that evolves in a manner 
consistent with nucleation-limited growth. Although the 
supersaturation is arbitrary and likely larger compared to 
the critical supersaturation in experiments, a systematic 
variation of the driving force allows us to extract the key 
features associated with the near-equilibrium behavior of 
the nanowire during its growth. The supersaturation for 
growth is generated by attempting to place Si atoms on 



the surface within an MD simulation and then accept- 
ing the move if an Au atom is within a nearest-neighbor 
distance. Alternate procedure wherein the droplet is su- 
persaturated by directly introducing atoms within an iso- 
lated droplet have negligible effect on the ensuing dynam- 
ics; the segregation at the surface and interface is fairly 
robust. The MD simulations employed are same as those 
for diffusivity calculations. 



Si(lll)-AuSu interface kinetics and step energetics 

The simulation cell contains of a crystalline Si(lll) 
substrate (8640 atoms) abutting a AuSi alloy (34560 
atoms). The dimensions of simulation cell is 14.14ao x 
14.69ao X 25.97ao, where ao is the lattice parameter 
at 873 K. The initial AuSi alloy with the desired con- 
centration is randomly structured and MD simulations 
are performed at 873 K. Periodic boundary conditions 
are applied along the in-plane directions with free sur- 
faces normal to the interface. The Si(lll) substrate is 
initially fixed and the AuSi alloy is relaxed for 0.1ns. 
The planar nucleation/melting simulations are carried 
by fixing the bottom bilayer of the substrate. The 
crystalline and liquid Si atoms are identified using the 
local order parameter described earlier. The simula- 
tion cell for extracting the step enthalpies is similar; 
it 19.80ao X 19.60ao x 10.39ao in size and consists of a 
Si(lll) slab (16128 atoms) and a AuSi alloy slab (16128 
atoms). The composition of the initial AuSi alloy is 
^Si = 0-46 and it is relaxed using isobaric, isothermal 
MD simulations (Supplementary Methods). The simu- 
lation cell is fully periodic so that the cell consists of 
Si( 111)- AuSi interfaces. The dimension normal to the 
free surface adjusts to maintain zero pressure while two 
in-plane dimensions (110) and (112) are fixed. Follow- 
ing relaxation (0.3ns), the total potential energy Ui is 
averaged over 0.1ns. Single bilayer steps are created on 
both surfaces of Si slab while preserving the atoms in 
the Si(lll) crystal and liquid AuSi. The steps are ori- 
ented along the index (112) direction. The total poten- 
tial potential energy is recorded as /72 and the difference 
(J = {U2 — Ui)/L = 0.16eV/nm, where L is the total step 
length. 

Si adatom energetics on truncating and main 
facets: Canonical MD simulations are performed at tem- 
perature T = 1 K. Each computational cell consists of 
an unreconstructed crystal Si slab, free at both both 
surfaces and periodic in the slab plane. The dimen- 
sions of simulation cells are [7.07 x 9.38ao x 6.63ao] and 
[7.07ao X 7.35ao x 5.20ao] for (113) and (111) surfaces and 
contains 3520 and 2160 atoms, respectively. The slabs 
were relaxed for 10 ps. A single Si atom is placed within 
the interaction range on each of the surfaces and relaxed 
for another 10 ps. All possible positions with maximum 
nearest neighbors were simulated to obtain the lowest po- 
tential energy. The maximum potential energy difference 
before and after addition of the Si adatom is recorded as 
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the Si adatom enthalpy Had- The extracted enthalpies 
are — 5.4eV/atom and — 3.5eV/atom for the (113) and 
(111) surfaces, respectively. Note that the stepped struc- 
ture of unreconstructed (113) surfaces can be interpreted 
as (111) steps separated by single atomic rows of (001) 
terraces. The {120} facets are also stepped and we there- 
fore expect similar trends. 
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FIG. 8: (top left) Atomic configuration and the corresponding normalized density profiles as a function of distance from the 
surface {h — H) in an AuSi thin film at the eutectic point for the model system, (a) The normalized density profiles at (a) 
eutectic point X^^(590K) = 0.31, (b) elevated temperature T = 873 K, and (c) elevated temperature and supersaturation, 
Xsz(873K) = 0.46. 




FIG. 9: a) Surface tension 7 versus — /3 for AEAM Au at 1400 K. /3 = —0.025 is a fit to the experimental value, (b) Segregation 
profiles for the AEAM Au/Si alloy thin film, before and after charge gradient correction. 
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FIG. 10: Temporal evolution of the ensemble- averaged mean square displacement of surface and bulk atoms in equilibrium MD 
simulation of an equilibrated 2R = 10 nm isolated Au53Si46 particle. 




FIG. 11: Equilibrium morphology of a 2R = 10 nm nanowire with {112} sidewall facets, (left inset) Schematic showing the 
relationship with 110 sidewall facets, (left) Atomic configuration of the complete nanowire-droplet-substrate system, (middle 
inset) (right) A mid-section slice along one of the (112) directions. Dashed lines in the figure indicate the truncating facets. 
The SiNW structure is the same as in the double-capped simulations, i.e. it has a hexagonal cross-section with {110} facets. 
The nanowire stem is 3.5 nm in length. 
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FIG. 12: Same as in Figure 5 in the main text, but for Xsi = 0.485 (At ~ 45 K). The t = 4ns configuration is shown with and 
without hquid Si (red) atoms. We see multiple growing nuclei on the main facet indicating that they are supercritical. The 
nuclei form preferentially near the corner edges where the main facet meets two truncating facets. 




FIG. 13: Atomic-configuration of an equilbrated eutectic droplet AussSie? on a Si(lll) surface. The color scheme is similar to 
that in the main text with the exception that we make no distinction between liquid and solid Si atoms. The inset shows the 
morphology of eutectic AuSi droplets on Si(lll) reported in prior experiments [P. The image is reproduced with the authors' 
permissions. 
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Supplementary Methods 
Atomic-scale simulations 

The semi-grand canonical Monte-Carlo simulations employed to relax the system at the growth temperature consist 
of translational and Au^Si exchange (transmutation) moves [2l[3]. The ratio between translational and exchange 
moves is fixed at 10 : 1. The acceptance ratio for translation moves is constrained to 50% by adjusting the magnitude 
of the translations. The equilibration is performed until the energy, pair distribution functions and density profiles 
converge. In some cases, canonical (NVT) Monte-Carlo computations are performed by turning off the transmutation 
moves. 

The MD simulations are performed by integrating the Newton's equations of motion using velocity- Verlet scheme. 
The time step is 1 fs and the trajectories for the canonical ensemble are generated using a Nose-Hoover to fix the 
temperature at 873 K [2 . Isobaric isothermal simulations (NPT) are performed by using a Parinello- Rahman barostat 
that allows the cell volume to adjust to the imposed press. The fictitious mass of the barostat is fixed at XX. 



Contact angle of 2R = 10 nm droplet on Si(lll) surface 

A 2R = 10 nm size spherical droplet consisting of ^ 32000 atoms is constructed by assigning random coordinates 
about positions corresponding to average Au-Si equilibrium distance based on the model angular EAM potential. 
One half of the droplet is then placed on an atomically flat Si(lll) substrate (thickness 8 bilayers). The substrate 
dimensions are as mentioned in the main text. It is periodic along the in-plane directions. The bottom layer is fixed to 
prevent the entire system from translating. The exposed regions are decorated with a {^/3 x a/3) 30° Au monolayer 
and then entire system is then equilibrated to the growth temperature using SGMC minimization followed by NVT 
MD. Upon equilibration, we see some evidence for disorder within the monolayer. Systematic thin film studies indicate 
that for our model system, the monolayer undergoes a surface transition at and around 900° K. A similar transition 
has been reported in experiments [2|, albeit at 1000°K. The difference between the experimental and computational 
contact angles is around 7°. We attribute this mainly to the differences in the manner in which the droplets attain 
the eutectic composition. In experiments, the silicon from the substrate melts into the initially Au-rich droplet to 
arrive at the equilibrium composition, resulting in a the crater like morphology of the substrate. On the other hand, 
the droplet in the computations is already at liquidus composition before it stabilizes on the substrate. The substrate 
is almost flat and therefore increases the measured contact angle with respect to the substrate surface plane. 



Charge gradient corrections for AEAM 

The computational cell for the MD simulations is chosen to be an FCC crystal slab with dimensions 10 ao x 10 ao 
in the periodic surface plane, and a thickness of 8ao. The slabs are melted and equilibrated for Ins; the effect of 
temperature was studied by equilibrating the system for 250 ps. Data runs 1 ns in duration are carried out for each 
temperature T. The mechanical definition is used to extract the surface tension [5], 

7 = 2 y [^^^ + "^yy - ^^^^1 ^ ^ Z^i + ^ ^ * ^ ^ 

-oo ^ ' 

Here the aa/s are components of the surface stress tensor, a, /3 = 1, 2, 3 along the x, y and z directions in a Cartesian 
coordinate system. S is the surface area, is the mass of atom Pa is the ath momentum component of atom z, 
is the ath component of the distance atoms i and j, fa is the ath component force component of atom and both 
the summations over i and j run over all the atoms in the systems. 

The Au-Au interactions employed as part of the AEAM framework used in this study are of the form suggested in 
Ref. [6]. For pure liquid Au at T = 1400 K, the extracted surface tension is 0.9 J/m compared to the experimental 
value of l.llJ/m. The difference is significant and can affect the segregation qualitatively. To test this, we employ 
a charge-gradient correction U to fit the Au-Au interactions to the experimental value of the liquid/vapor surface 
tension; the Au-Si and the Si-Si interactions remain unchanged. 
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In this framework, the potential energy Ui of atom i is, 



rj,=pi+l3\Vpi\\ pi = J2nr). (S2) 

Here F is the usual embedding potential energy, ^ is the pair potential energy, pi the local charge density at atom z, 
and r]i is the local charge density of atom i with gradient corrections, and r is the distance between atoms i and j. 
The constant (3 emerges as the main fitting parameter for the liquid surface tension. Notice that the formulation in 
Ref. [7] introduces another fitting parameter c which is set to zero here. Then, the ath component of force rriiaf on 
atom i becomes, 

=-2:[2d^ + ^ V (S3) 

* \3i^^ 3i^^ j¥« j¥« 

where the first term is the standard EAM-based force and the second term is the charge gradient correction contri- 
bution. Using this formulation, we have re-calculated the surface tension of pure liquid Au 7 as a function of the 
fitting parameter /3 at 1400 K. The results are shown in Supplementary Figure S2. The surface tension increases 
monotonically with the magnitude of /3 with the value /3 = —0.025 in excellent agreement with the experimental 
value. Simple tests are first performed to quantify the effect of this modified potential on the bulk properties of Au 
at 1400 K using MD simulations. For /3 = —0.025, the cohesive energy is — 3.626 eV/atom and the lattice constant is 
ao = 4.197^4. As comparison, for /3 = which reduces to the original AEAM potential, the corresponding values are 
— 3.624 eV/atom and ao = 4.206 A, respectively. 

This modified Au-Au potential within the AEAM formulation is employed to recalculate the density profiles for 
AuSi thin films. The surface segregation in AuSi film with 33% Si for charge gradient correction /3 = 0.025, which 
is also qualitatively same the result without charge gradient correction. The ensemble-averaged segregation profiles 
are show in Fig. |9^. Again, we see strong subsurface Si segregation which is decorated by a submonolayer of Au. 
The quantitative differences are minimal, notably that the peak associated with Au submonolayer decoration on the 
surface is slightly decreased. 



Supplementary Discussion 1 

We estimate the nucleation barrier on the planar Si(lll)-AuSi interface by assuming that the nuclei is a circular 2D 
island. The barrier for homogeneous nucleation on the planar Si(lll) surface is roughly AG ~ a'^h/ {p jisi)^ where <j 
is the step energy per step step (bilayer) height /i, p is the density of solid silicon, and psi is the free energy difference 
between solid and liquid phases. Our extracted values of the step enthalpies on Si(lll) vary between ah ^ 0.16 eV/nm. 
jisi for the small supersaturations in the AuSi model system is ~ lOOA^kJ/mol, approximated as the enthalpy of 
mixing predicted by the model AuSi system [34]. Using these values yields a barrier AG ~ 20 me V/ Ax that is of the 
same order as the extracted value of AG ^ 7meV/Ax. The self-consistent comparison is strong evidence that the 
nucleation barrier during nanowire growth is lowered due to the multi- faceted morphology of the growth front. 



Supplementary Discussion 2 

The Si incorporation rate follows from the volume growth rate due to the flux made available by the surface catal- 
ysis. It scales as kyipAg^ where p is the precursor gas partial pressure, kyi is normalized catalyst rate constant per 
unit partial pressure, and As is the exposed surface area of the droplet. The rate constant reported in recent experi- 
ments on growth of Si nanocrystals within an AuSi droplet exposed to a silane flux is kyi ~ 10 nm s -^Torr"^ [26]. 
Approximating the droplet surface area to that of a hemisphere of radius 5 nm, under typical processing conditions 
{p ^ 10~^ Torr) yields the incorporation rate assumed in the text. As a simple check, the axial nanowire growth rate 
then is of the order of a nanometer per minute which is consistent with growth rate observed in in-situ experiments. 
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The transition from flux-controlled to nucleation-controUed growth occurs when the flux due to catalysis equals the 
crystallization rate on the main facet, 

From Fig 6b in the main text, we have vq = 0.58 m/s and AG/ksT = 0.055. For Si volume flux of the order of 
I Si ~ 1 atom/ms, the axial nanowire growth rate is, 

- (S5) 



where p = 50 atom/nm^ is the density of crystallized Si and R = 5nm is the nanowire radius in our simulation. 
Equating Eqs. SI and S2, we get = 0.002. 
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